Effect of phase relaxation on quantum superpositions in complex col·lisions 
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We study the effect of phase relaxation on coherent superpositions of rotating clockwise and 
anticlockwise wave packets in the regime of strongly overlapping resonances of the intermediate 
complex. Such highly excited deformed complexes may be created in binary collisions of heavy 
ions, molècules and atòmic clusters. It is shown that phase relaxation leads to a reduction of the 
interference fringes, thus mimicking the effect of decoherence. This reduction is crucial for the 
determination of the phase-relaxation width from the data on the excitation function oscillations in 
heavy-ion collisions and bimolecular chemical reactions. The difference between the effects of phase 
relaxation and decoherence is discussed. 
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I. INTRODUCTION 

For highly excited strongly interacting many-body sys- 
tems the independent particle picture has very limited 
validity 0. For high excitations, the interaction results 
in a quick decay of single-particle 0, EI H as well as 
collective modes 0, 0] , which are not eigenstates of the 
Hamiltonian of the system. This decay leads to the for- 
mation of complicated many-body configurations. Each 
of these many-body states is characterized by a uni- 
form occupation of the accessible phase space and shar- 
ing the energy between many partides of the system. 
The characteristic time for the formation of such spa- 
tially extended ergodic many-body states, independent 
of the initial conditions, is given by the inverse spreading 
width, T crg = h/T spl - 4]. Introduced by Wigner [jj, T spr 
also characterizes the width of the distribution of the ex- 
pansion coefhcients of the many-body eigenstates over a 
noninteracting mean-field basis. For sufficiently high ex- 
citation energy the decay of single-particle modes creates 
a linear superposition of a very large number of many- 
body configurations. The question is whether the phase 
relations between these individually ergodic many-body 
states in the superposition may still preserve a mem- 
ory of the way the system was excited. This question 
is of fundamental importance for the study of relaxation 
phenomena in nuclear, atòmic, molecular and mesoscopic 
many-body systems, and for many-qubit quantum com- 
putation. In particular, if phase relaxation is longer than 
energy relaxation h/T'apr, this can extend the time for 
quantum computing |g, Q beyond the quantum-chaos 
border [10. 

To answer this question from first principies is difhcult 
due to computational limitations. Indeed, in order to 
solve the full quantum many-body problem one would 
rcquirc a many-qubit quantum computer. Therefore, 
the only currently available resort to search for possi- 
ble manifestations of long phase relaxation is the exper- 
iment, and a careful data analysis. Nuclear systems are 



an ideal laboratory to study many-body systems, since 
nuclear interactions are so strong that external pertur- 
bations can often be neglected. In particular, the anal- 
ysis of Refs. [Ü 0, E3 of the data on asymmetry around 
90° cm. in angular distributions of evaporated protons 
from heavy nuclei in nucleon-induced and photonuclear 
reactions clearly indicates that phase relaxation time in 
highly excited strongly interacting many-body systems 
can be up to eight orders of magnitude 7] longer than 
energy relaxation time h/T spr . This reveals a new form of 
matter-thermalized non-equilibrated matter introduced 
by one of us in Refs. [Ulli El. 

A more subtle indication of a slow phase relaxation is 
found in light heavy-ion scattering [3> El Eï E3 • It 
manifests itself in the formation of rotating wave pack- 
ets whose spreading time, which is given by the inverse 
phase-relaxation width, is much longer than the energy 
relaxation time 7ï/r S p r . The approach 0,0,0] treats 
time-delayed light heavy-ion scattering in terms of the 
formation and decay of quasimolecular resonances |l8j . 
The highly excited coherently rotating intermediate sys- 
tem has the energy of the intrinsic excitation > 10 MeV. 
The intermediate system is in the regime of strongly over- 
lapping resonances. Therefore, this coherent rotation 
does not originate from the discreteness of the spectrum, 
which is not resolved, but is due to the spin off-diagonal 
correlations between partial width amplitudes. Indeed, 
the period of the coherent rotation is much shorter than 
the inverse level spacing of the intermediate system. This 
reveals a new root, as compared to Bohr's correspondence 
principle, for a quantum-classical transition in highly ex- 
cited many-body systems. 

The width of the rotating wave packets is about 
í/d + f3t/h w l/d + /36/huj, where d is the effective num- 
ber of partial waves, (3 is the spin off-diagonal phase- 
relaxation width, lü is the angular velocity, and t and 
are the time and angle of the rotation, respectively. 
This seemingly allows us to determine the wave-packet 
spreading rate from the time power spectra at differ- 
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ent scattering angles, which can be reconstructed, for 
binary collisions, from the data on energy fluctuations 
of the cross sections 0, 0|. These fluctuations origi- 
nate from the energy-fluctuating collision amplitude cor- 
responding to a resonance time-delayed process. How- 
ever, for elàstic heavy-ion scattering, such energy fluctu- 
ations can be typically observed only at backward angles 
140° < 8 < 180°, since for smaller 8 the direct reac- 
tion contribution grows exponentially and becomes much 
greater than the time delayed one. Therefore, since at the 
initial moment of time the deformed intermediate system 
is oriented in the forward direction 8 ~ 0°, the spreading 
of the wave packets within the backward angular interval 
140° < 8 < 180° can hardly be detected reliably. More- 
over, the width of the wave packets at this backward- 
angle range does not allow to determine (3, since the an- 
gular dispersion 1/d at t = is also unknown. In this pa- 
per we show that both d and [3 can indeed be dctcrmined 
unambiguously. This results from the strong sensitivity 
of the interference fringes between the rotating clockwisc 
and anticlockwise wave packets for the backward angle 
range to the the phase relaxation width: As (3 grows, the 
interference fringes are suppressed more strongly. The 
effect is considered in relation to heavy-ion scattering 



and bimolecular chemical reactions. The difference be- 
tween the effects of phase relaxation and decoherence is 
discussed. 



II. THE TIME POWER SPECTRUM OF THE 
TIME-DELAYED COL·LISIONS 

Following Ref . > we consider spinless collision part- 
ners in the entrance and exit channels. The time and 
angle dependent intensity of the decay (the time power 
spectrum), P(t,8), is given by the modulus square of 
the Fourier component of the energy-fluctuating colli- 
sion amplitude. It can be also expressed as the Fourier 
component of the amplitude energy autocorrelation func- 
tion. P(t, 8) has been obtained in |l9j | by summing over 
a very large number of strongly overlapping resonance 
levels, T/D 3> 1, where T is the total decay width and 
D is the average level spacing of the intermediate com- 
plex. As a result, after changing from summation over 
the resonance levels to integration, which is a good ap- 
proximation for t <C h/D, P(t, 0) takes the form: 



P(t, 8) oc H(t) exp(-rt/n) Y,( 2J+ !)( 2J ' + 1)[W(J)W(J')} 1/2 exp[i($-ujt)(J-J')-/3\J-J'\t/n)Pj(d)Pj>(6). (1) 
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Here, H(t) is the Heaviside step function, fi is the 
spin phase-relaxation width, u) is the angular velocity 
of the coherent rotation, $ is the deflection given by 
the total spin J (in h units) derivative of the poten- 
tial phase shifts (direct reaction), and Pj(0) are the 
Legendre polynomials. The physical meaning of the in- 
verse spin phase-relaxation width, %/ f3, is the character- 
istic time for the angular spreading of the clockwise and 
anticlockwise rotating wave packets. The partial aver- 
age reaction probability is taken in the J-window form, 
W(J) = {\SS J (E)\ 2 ) oc exp[-(J - J) 2 /d 2 }, where J is 
the average spin and d is the J-window width. 

It should be noted that a similar effect of coherent rota- 
tion of the nuclear molecule, as described by Eq. (JIJ with 
fi = 0, was found in Ref. [2(| for a scattering of heavy nu- 
clei. However, while Eq. Q has been obtained for a large 
number of open channels and strongly overlapping reso- 
nances of the intermediate system |l9j, the time power 
spectrum in Ref. |20j was derived from one open channel 
and one resonance pole form of the S'-matrix elements. 
Therefore the latter approach is not applicable for the 
fortheoming analysis of the 12 C+ 24 Mg system since, for 
this system, the number of open channels is much greater 
than one. Furthcrmore, the one-resonanec pole form of 
the S'-matrix elements employed in Ref. |20f results in iso- 



lated resonances with the energy spacing between them 
to be about hus. This is not the case for the 12 C+ 24 Mg 
elàstic scattering analyzed in this paper. Indeed, the 
data on the excitation function for the 12 C-(- 24 Mg elàs- 
tic scattering at 8 — 180° reveal about 15 local màx- 
ima on the energy interval AE cm = 13 — 22 MeV pi) . 
Then the interpretation [20| would mean that average 
level spacing, ~ 0.6 MeV in our case, between the lo- 
cal màxima is given by Tiu. This, in turn, would re- 
flect the quasiperiodic behaviour of the cross section en- 
ergy autocorrelation functions with the period of about 
0.6 MeV. This is inconsistent with the data [2Ï| showing 
the quasiperiodicity of the cross section with the period of 
about 2.9 MeV yielding hu = 1.35 - 1.45 MeV 0113. 
Moreover, the interpretation [2(| implies that these lo- 
cal màxima in the excitation function should be approx- 
imately equidistant. Instead, the data show big 

fluctuations, from 0.3 MeV to 1.3, for the energy spac- 
ing between the nearest-neighboring local màxima in the 
excitation function. Note that the indication against 
the interpretation in terms of isolated resonances comes 
from the statistically insignificant channel-channel corre- 
lations for the 12 C+ 24 Mg elàstic and inelastic scattering 
at 8 = 180° 0| . Finally, the effect of coherent rota- 
tion does not survive a generalization j23 by taking 
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FIG. 1: (Color online) Angular dependence of the decay 
intensity of the highly excited intermediate complex at dif- 
ferent moments of time for different combinations of P and 
d. Panel (a) corresponds to time í = 3T/8; (b) t = 7T/16; 
(c) í = T/2; (d) t = 5T/8. Here, T is the period of one com- 
plete revolution of the complex. Sòlid lines (black) correspond 
to = 0.01 MeV and d = 3; dashed (red) to p = 0.04 MeV 
and d = 4; dashed-dotted (blue) to = 0.075 MeV and d = 8. 
The bigger /3 the smaller the interference fringe contrast is. 



into account intermediate structure by means of inclusion 
of many isolated resonances in the 5-matrix elements. 
On the basis of the above arguments we believe that the 
interpretation of the fine structure in the excitation func- 
tions for the 12 C+ 24 Mg scattering at 9 = 180° in terms 
of overlapping, rather than isolated, resonances of the 
intermediate system is a justified approach. 

First, we calculate P(t, 9) for the set of parameters 
obtained from the description 17] of the experimen- 
tal cross section energy autocorrelation functions [2Ï| 
for 12 C+ 24 Mg elàstic and inelastic scattering at 9 = 
180° . For these col·lisions the analysis of the oscil- 
lations in the cross section energy autocorrelation func- 
tions 0, 0| indicates the formation of stable rotational 
wave packets, in spite of the strong overlap of resonancc 
levels in the highly excited intermediate molecule. The 
set of parameters is [ï^: $ = 0, d = 3, J = 14, 
P = 0.01 MeV, Hlü =1.45 MeV, and T = 0.3 MeV. 




FIG. 2: (Color online) Time dependence of the decay inten- 
sity of the highly excited intermediate complex at (a) 9 = 180° 
and (b) 9 — 170.6°, for the different sets of and d. Sòlid 
lines (black) correspond to P = 0.01 MeV and d = 3; dashed 
(red) to = 0.04 MeV and d = 4; dashed-dotted (blue) to 
P = 0.075 MeV and d = 8. The smaller p the deeper the min- 
imum is, due to the interference between the wave packets for 
9 = 170.6° at t = T/2. 



In Figs. U] we plot the quantity AP(t,6)/(a(E,6)) at 
four moments of time, with T = 2n /u being the period of 
one complete revolution. Here, (a(E,6)) oc J °° dtP(t,ff) 
is the energy-averaged differential cross section for the 
time-delayed collision. In Figs. ^ P(t, Q) is scaled with 
(cr(E, 9)) for the reasons discussed in Ref. The con- 
stant A is determined by the condition AP[t = 0, 9 = 
Q)/{a(E,9 = 0)) = 1. Figure[H» shows the two slightly 
overlapping wave packets rotating towards each other in 
the backward direction. In panels (b) and (c), the wave 
packets strongly overlap around 9 = 180°, producing in- 
terference fringes. Finally, in panel (d) the wave packets 
have passed each other and move apart rotating in the 
forward direction. 

In Figs. UI we also plot AP(t,9)/{a(E,9)) with the 
same constant A and the set of parameters as before, 
except for d and 0. We consider two other cases: d = 4, 
p = 0.04 MeV, and d = 8, P = 0.075 MeV. One can see 
that when the wave packets overlap only slightly (pan- 
els (a) and (d)) one can hardly distinguish between the 
three combinations of d and j3, even though /3 changes 
by almost one order of magnitude. On the other hand, 
from panels (b) and (c), we observe that a fringe con- 
trast due to the interference of the wave packets is very 
sensitive to absolute vàlues of P and d. This result is fur- 
ther illustrated in Fig. 2 where we plot P(t, 9) calculated 
for the three sets of j3 and d vàlues, as a function of t for 
9 = 170.6°. One observes that the smaller P is the deeper 
minima are, due to the destructive interference between 
the wave packets at t = T/2. However, for 9 — 180°, 
P(t,9) in Fig. |3 is insensitive to the three sets of P and 
d. 
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Note that in order to reproduce the width of the wave 
packets in Figs. ^ and [21 for the maximal possible d = I 
value, one has to take /3 = 0.1 MeV. Therefore, even the 
largest possible [3 value is more than one order of mag- 
nitude smaller than r spr . The latter is usually estimated 
from the width of giant resonances 0, [1 : r spr ~ 5 MeV. 
This indicates that energy relaxation, i.e., the process of 
formation of ergodic many-body configurations, in highly 
excited nuclear systems can be much faster than phase 
relaxation. 



III. DETERMINATION OF THE TIME POWER 
SPECTRUM FROM THE DATA ON THE 
EXCITATION FUNCTION OSCILLATIONS 

We consider a collision of the spinless reaction partners 
in the entrance and exit channcls. The cross section of 



is the instantaneous decay amplitude received at a de- 
tector at time t with the time uncertainty of about h/I, 
for a collision initiated at t = 0. In the sum in Eq. J2Jl, 
N = I/AE and (N + 1) is the number of equidistant en- 
ergy vàlues for which the detailed angular distributions 
should be measured. The summation in the r.h.s. of 
Eq. J2J may be changed to the integration if we use a lin- 
ear interpolation for f(E, 9) in between the conseqüent 
energies E n and E n+ i. This is justified provided the en- 
ergy step AE is smaller or about the minimal character- 
istic energy scale of a variation of the cross section. For 
example, in Ref. [2^ the energy step is AE cm = 133 keV 
while the minimal characteristic scale of the variation of 
a(E,9= 180°) is T = 0.3 MeV jHG]!- 

Having reconstructed V(t, 9) from the excitation func- 
tion data one obtains the time power spectra P(t, 9) = 
\V(t,9)\ 2 for any angle. The later quantity is equiva- 
lent to that directly measured using real-time methods 
of femtochemistry, with the energy bands of / for the 
pump and probe làser pulses, to monitor chemical re- 
actions dynamics |2^,[2g. Due to its dependence on 9, 
P(t, 9) determined by the proposed method could provide 
even more detailed information on the collision dynam- 
ics, as compared with the femtochemistry experiments, 
since latter measurements often allow to obtain P(t, 9) 
averaged over 9. 

For the energy interval / being long enough to pro- 
vide a sufRcient time resolution h/I to resolve interfer- 
ence fringes [2?} the proposed experiment would allow 
to test our theoretical predictions. If these were con- 
firmed by the data, then the experimentally determined 



the collision is given by a(E,9) = \ f(E : 9)\ 2 . Here E is 
the energy in the entrance channel and 9 is the scattering 
angle. The reaction amplitude f(E, 9) is a linear combi- 
nation of S'-matrix elements S J (E) corresponding to the 
different total spin vàlues J. Suppose that, for a fixed 
energy E\ 1 a measurement of detailed angular distribu- 
tion is performed. Then, using, e.g., the method j^, one 
can find S J (Ei) and, therefore, f{E\,ff). Measuring the 
angular distributions for different energies on the energy 
interval / = E min — E max with the energy step AE and 
energy resolution smaller or about the minimal charac- 
teristic energy scale of a variation of f(E, 9) one obtains 
a detailed energy and angle dependències of the collision 
amplitude. Its Fourier component, 



(2) 



P(t, 9) can be compared with our calculations presented 
m Figs. □ and El with the purpose of evaluation [3 and d. 
In particular, the energy interval of 10.5 MeV, over which 
excitation function for the 12 C+ 24 Mg elàstic scattering 
at 6 = 180° was measured [%|, should be sufRcient to 
resolve the interference fringes in P{t, 9) |2Ï| . 

It should be noted that the method proposed above al- 
lows to reconstruct P(t, 9) which includes both potential 
scattering, i.e. fast processes, and the time delayed mech- 
anism. Since the previous analysis [Ï6l IÏ7l| demonstrated 
a presence of a strong potential scattering component 
we expect that this method will yield a strong potential 
scattering contribution in P(í, 9). However this contribu- 
tion is expected to be restricted to relatively short times 
tp S < h/AE ps , where AE ps is a characteristic energy in- 
terval for the variation of the potential scattering ampli- 
tude. For AE ps ^ 2-3 MeV [l·|, t ps ~ 3-2 x 10" 22 sec. 
Therefore, this potential scattering contribution is ex- 
pected to be restricted to t/T < 0.1 in Fig.|2and should 
not show up at all in Fig. ^ 

The above method to obtain P(t, 9) from the data 
is model independent and, therefore, most reliable one. 
However it is very demanding experimentally since 
a(E, 9) must be measured for a very large number of 
angles 24]. There is another way to extract P(t, 9) from 
the data even if the excitation function is measured for 
a single angle. This method can be applied if the di- 
rect reaction process provides a major contribution to 
the cross section 0. We consider again a case of spin- 
less collision partners in the entrance and exit channcls. 
The main idea is based on the decomposition of the col- 
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V(t, 9) oc / dEexp(-iEt/h)f(E, 9) w AE V exp[-i(£ min + nAE)t/ti]f{E mm + nAE, 9) 



71=0 

r 
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lision amplitude, f(E, 9) = fdir(E, 9) + Sf(E,9). Here 
fdir{E, 9) is the energy smooth amplitude corresponding 
to potential scattering or direct fast process, and Sf(E, 9) 
is the energy fluctuating amplitude, (6f(E,9)) = 0, 
corresponding to the time-delayed collision. Considcr 
first an idcalized case of energy independent amplitude 
fdir{E, 9) when it coincides with the energy averaged am- 
plitude (f(E,9)}. The collision cross section has the 
form a(E,9) = \f(E,9)\ 2 = a dir (E,9) + a fl (E,9) + 
2Re[6f(E,9)f dir (E,9)*]. Here the potential scattering 
or direct reaction cross section a dir (9) = \fdi r (E,9)\ 2 is 
energy-independent and <jfi(E,9) = \Sf(E,9)\ 2 . For a 



relative contribution of the potential scattering or direct 
reaction cross section to be about 70% or more of the 
total cross section, Aa(E,9) = a(E,9) - (cr(E,9)) ~ 
2Re[Sf(E,9)f dir (E,9)*}, where (a(E,9)) = a dir (9) + 
(afi(E,9)} is the energy averaged cross section. This 
means that ràpid energy variations of the cross section 
originate mostly from interference between the energy 
fluctuating 5f(E, 9) and the energy independent ampli- 
tude f dir (E,9). Suppose that the excitation function is 
measured on energy interval / = E max — E m i n with en- 
ergy step AE so that the total number of steps is N + 1 , 
where N = I/AE. Then, for t > 0, we have [13 



dEexp(-iEt/h)Aa{E,9) oc ( dE exp{-iEt/h)Sf(E,9) 

'min J E min 

N 



AE J2 cxp[-i(E mm + nAE)t/h]Sf(E min + nAE, 9) <x V(t, 9) (3) 
I 



n=0 



with P{t,9) = \V(t,9)\ 2 . In Eq. Q we have employed a 
causality condition, P(t < 0, 9) = 0. This condition im- 
plies that the molecule cannot decay before it is formed 
at t — 0. Again the above information on P(t, 9) is equiv- 
alent to that obtained directly using real-time methods 
of femtochemistry to monitor unimolecular chemical re- 
actions, with the energy bands of / = {E max — E m i n ) for 
the pump and probe làser pulses [2í| . 

The above consideration can be extended to a case 
when the potential scattering or direct reaction ampli- 
tude fdir{E,9) depends on energy but this dependenec 
is considerably smoother than the energy dependenec 
of the amplitude Sf(E, 9) for the time-delayed colli- 
sion. We shall again assume that a relative contribu- 
tion of direct processes into the energy averaged cross 
section is about 70% or more. We define the charac- 
teristic energy interval Id of a variation of fdir (-E 1 , 9) 
with íj > T, where T is the characteristic energy in- 
terval of variation of 5f(E,9). The energy interval 
I d can be evaluated using, e.g., the trend reduction 
method [!|. In particular, for the 12 C-(- 24 Mg scat- 



tering 22] analyzed in Sec. II, this interval was eval- 
uated to be about 4 — 5 MeV |2Ï|. i.e. much greater 
than r = 0.3 MeV. Then one can find the best poly- 
nomial fit of the energy dependence of the cross section 
a{E, 9) = c7 d (E, 9) + a fl (E, 9) + 2Re[Sf(E, 9)f dir {E, 9)*} 
with the order of polynomial being [I/I d ] + 1, where [I/I d ] 
is integer part of I / I d The resulting energy smooth cross 
section is denoted as à(E,9) ~ a d (E,9) + (afi(E,9)} 
with {(jfi(E,9)) being energy averaged, i.e. energy in- 
dependent, cross section of the time-delayed processes. 
It can be found using the Standard procedure em- 
ployed for the analysis of Ericson fluctuations (see, e.g., 
Refs. j^H^I). As a result, one finds energy smooth di- 
rect reaction cross section ad(E,9). Consider a(E,9) — 
à(E, 9) ~ 2Re[Sf(E, 6)f dir (E, 9)*}. Assume first that en- 
ergy dependence of the phase of fdir(E,9) is negligible, 
arg[fdi r (E,9)] — constant, and the smooth energy de- 
pendence of <Jd(E, 9) is due to the energy dependence of 
\fdir(E,9)\ = a d (E,9) 1 / 2 . Then applying the same ar- 
guments as those used to obtain Eq. for t > 0, we 
have 



J 



JV 

P(t, 9)(xJ2 exp[-i{E min + nAE)t/h] [a(E mm + nAE, 9) - a(E min + nAE, 9)}/a d (E min + nAE, 9) 1 ' 2 (4) 

71=0 



with P{t,9) = \P{t,9)\ 2 . 

In order to generalize Eq. by taking into ac- 
count energy dependence of arg[f d i r {E, 9)] = (f>(E,9) 
(f dir (E,9) = \f dir (E,9)\exp[i<t>{E,9)]) we use the 



linear approximation <fi(E, 9) — 4>{E, 9) + (E — 
É)d(j}{E,9)/dE\ E=È1 where É = (E max - E min )/2 and 
d(j)(E,9)/dE\ E= E = tdirj'h with tdi r <C h/T being the 
time delay due to the potential scattering or direct re- 
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actions |29j. Then, it is easy to see, that, for t > tdi r , 
the r.h.s. of Eq. is changed to V(t - t dir ,9). This 
manifestly demonstrates that the interference between 
time-delayed and fast direct processes in the cross sec- 
tion is a precondition for monitoring the time evolution, 
the fast process switches on the clock at the initial mo- 
ment of time td playing the role of the pump pulse. In 
the absence of direct processes the initial moment of time 
is not defined. 

Note that P(t, 9) can also be expressed as a co- 
sine half-Fourier transform of the cross-section cn- 
ergy autocorrelation function C(e,9), provided the rel- 
ative contribution of potential scattering is larger than 
70% [HG3- In Fig.Elwe plot C(e,9)/C(e = 0,9) = 
[J™ dtcos(et/h)P(t,9)}/(a(E,9)) for (a) 9 = 180° and 
(b) 9 — 170.6°, for the three combinations of (3 and d. 
Since the P(t, 8 = 180°) are close for the different sets 
of (3 and d, the corresponding C(e, 9 = 180°) can hardly 
be distinguished reliably. This is the reason that the 
analysis at 6 = 180° pj| E3 is not sufficient to deter- 
mine unambiguously the vàlues of (3 and d. However, 
C(e,9 = 170.6°) is more sensitive to the different sets of 
(3 and d. In particular, for e > 4 — 5 MeV, the oscilla- 
tions in C(e,9 = 170.6°) for (3 = 0.01 MeV, d = 3 and 
f3 = 0.075 MeV, d — 8 are out of phase with the absolutc 
value of their diffcrcncc being up to 0.4. Still, from the 
comparison of Figs.^c),[2ïb) and Fig.[3]we observe that 
the sensitivity of P(t, 9) to different vàlues of (3 and d 
is considerably stronger than the sensitivity of C(e, 9) to 
these vàlues. Also P(t, 9) obtained using Eqs. ® and Q 
does not acquire statistical errors due to the finite energy 
interval /. Indeed such a method provides information 
on P(t, 9) equivalent to that obtained in the femtochem- 
istry experiments. Thcrcforc the only uncertainty of this 
method is a finite time resolution, Aí ~ fi/I, which is 
the uncertainty in the femtochemistry experiments with 
the finite energy band / of the pump and probe làser 
pulses jlMEBl- On the other hand, a reconstruction of 
P(t, 9) from C(e, 9) may result in additional errors due 
to possible statistical uncertainties in C{e,9) related to 
the finite energy range / 30] . Such additional statistical 
uncertainties may not be excluded even though the ef- 
fects studied in this paper can clearly not be associated 
with Ericson fluctuations. Indeed, unlike, e.g., oscillating 
structures in Fig. [3J Ericson theory prediets angle inde- 
pendent Lorentzian shapes for the C(e, 9). Therefore, 
for the reasons described above, we suggest that a recon- 
struction of P(t, 9) directly from the excitation functions 
(Eqs. Q and l@J) is more reliable than from C(e,9). 



IV. FROM HEAVY-ION COL·LISIONS TO 
BIMOLECULAR CHEMICAL REACTIONS 

An important question is whether the rotating wave 
packets and their interference can oceur in colliding sys- 
tems for which more rigorous approaches can be im- 
plemented. The encouraging answer has been given by 
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FIG. 3: (Color online) Directly measurable cross section en- 
ergy autocorrelation functions calculated at (a) 9 = 180° and 
(b) 9 — 170.6° for the different sets of f3 and d. Sòlid lines 
(black) corresponds to f3 = 0.01 MeV and d = 3; dashed 
(red) to f3 = 0.04 MeV and d = 4; dashed-dotted (blue) to 
13 = 0.075 MeV and d = 8. 



the calculations of the time evolution of the H + D2 — > 
H D + D bimolecular chemical reaction [3ÏI 133 . l33| . The 
evidence for rotational wave packets has been demon- 
strated (see, e.g., Fig. 3 in l31| ) due to the interference of 
the overlapping resonances j34j of the quasibound com- 
plex with J = 15 — 20 total spin vàlues, (i.e., within 
our approach, d ~ 5). In this peripheral time-delayed 
collision the rotational wave packets are stable show- 
ing a moderate spreading during half a period of the 
rotation. Within our description, this reveals a small 
(3 value in the time-delayed peripheral reaction mech- 
anism, i.e. Tr/3/huj 4; 1/d. Then, in accordance with 
our prediction, the coherent superposition of the rotating 
wave packets should produce a strong interference con- 
trast. This is indeed the case as clearly seen in Fig. 3 of 
Ref. [33 • Similar effects of the stable rotating wave pack- 
ets and their coherent superpositions with strong inter- 
ference contrast have also been found in the calculations 
for F + HD -> HF + D bimolecular reaction [HIH. 



The above results should allow to test our suggestion 
for the determination of the phase relaxation time and, 
in general, for the time evolution of the collision in bi- 
molecular chemical reactions. For example, in 32] a 
total time dependent amplitude was calculated for the 
H + D2 — > HD + D reaction. Then it was splitted 
into short-time direct reaction component and the time- 
delayed one. Taking the half-Fourier transform, explicit 
energy and angle dependències of both the direct and 
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time-delayed amplitudes were obtained. It was found 
that in the region of backward angles 9 < 70° and for 
the cm. energy interval 1.6-2 eV the direct reaction 
cross section is considerably bigger than the time-delayed 
one. Then one can apply our method (Sect. III) to ob- 
tain P(t, 9) of the time-delayed reaction mechanism from 
the energy dependence of the double-differential cross 
sections. This extracted P(t, 9) can be compared with 
its actual exact form since the latter is known and was 
used to generate the energy dependence of the double- 
differential cross sections. 

The calculations [U H3 suggests that in the forward 
angle range, where the Schròdinger cat states originated 
from the interference of the rotating towards each other 
wave packets are observed, contribution of the direct 
reactions is negligible. Therefore, for this particular 
H + Z?2 — > HD + D reaction for the cm. energy less 
than 2 eV our method can not be applied. Yet one can 
still check the sensitivity of the method by artificially 
adding an energy smooth background amplitude to the 
calculated one for the time-delayed reaction mechanism. 
The energy smooth amplitude would play a role of the 
direct reaction amplitude. This would enable one to di- 
rectly test the accuracy of the determination of the in- 
terference fringe contrast from the data on the double 
differential cross sections generated with the artificially 
added energy smooth direct reaction amplitude. Such a 
test may be important since for the higher energy and/or 
for some other bimolecular chemical reactions, a strong 
contribution of direct reactions is possible for the same 
angle range where the Schròdinger cat states show up. 

Note that pure energy resolution measurements of the 
cross sections of bimolecular chemical reactions, in par- 
ticular H + D 2 -> D + HD 37] and F + HD -> 
HF + D Q reactions, have recently become possible. 
Our results suggest that this should enable experimen- 
talists to extract information on the time dependence, in 
particular on the rotational coherent dynamics and possi- 
ble wave packet interference, of the bimolecular chemical 
reactions which previously could only be studied by the 
pump probe làser pulses technique to monitor unimolec- 
ular chemical reactions [2^, [2(j . 



V. PHASE RELAXATION VERSUS 
DECOHERENCE 

Since phase relaxation results in washing out the inter- 
ference fringes and eventually, for a sufficiently large /?, 
can destroy the Schròdinger cat states of Figs. ^ and El 
this effect appears to be quite similar to decoherence [3!j . 
Yet, there is an essential difference between the two ef- 
fect s. It is instructive to compare the extreme case of very 
fast phase relaxation (very large (3) with very fast deco- 
herence. Very fast phase relaxation corresponds to the 
regime of random matrix theory, i.e., to random phase re- 
lations between all partial width amplitudes correspond- 
ing to different strongly overlapping resonance states |^ . 



This yields angle independent exponential form for the 
time power spectrum P(t, 9) and unity for the normalized 
variance of the fluctuating cross section. Analogous to 
the fast phase relaxation would be quick decoherence be- 
tween all many-body quasi-stationary resonance states, 
similar to the decoherence of Fock states due to the cou- 
pling of the system to a phase reservoir [3íJ . Such a deco- 
herence would destroy interference terms between all the 
many-body eigenstates, leading to the vanishing of in- 
terference between different partial width amplitudes. It 
can be shown that in this case the time power spectrum 
would have the same angle independent exponential form 
as for fast phase relaxation. However, in contrast to fast 
relaxation, quick decoherence would result in washing out 
the cross section energy fiuctuations, reducing its normal- 
ized variance to D/r < 1 |40j . It should be noted that 
the suppression of electron transmission intensity fiuctu- 
ations through nanostructures due to dephasing, an effect 
similar to decoherence, can be described on the basis of 
different arguments allowing to evaluate dephasing rates 
from the amplitude of conductance fiuctuations j4Ï| . 



VI. DISCUSSION AND CONCLUSIONS 

Our approach is based on a new idea of slow phase re- 
laxation between partial width amplitudes in a regime of 
strongly overlapping resonances. Contrary to the conven- 
tional idea of random phase relations in this regime 0, 
IÍ2l |44| , our analysis invokes spin off-diagonal phase cor- 
relations and their slow decay for the resonance time- 
delayed col·lisions. Though these effects may oceur for 
atòmic cluster col·lisions [Ï3 and are already revealed for 
the bimolecular chemical reactions [^ÏLI^I^&lsEI^, the 
present consideration is essentially motivated by the pure 
energy resolution data on fiuctuations in 12 C+ 24 Mg elàs- 
tic and inelastic scattering at 9 = 180° [2Ïll22l|. demon- 
strating both the fine energy structure of ~ 0.3 MeV 
and quasiperiodicity with period ~ 3 MeV [ïil IÏ7| . It 
is these peculiar features that led us to the interpreta- 
tion of the 12 C+ 24 Mg scattering in terms of highly ex- 
cited quasimolecular states with strongly overlapping res- 
onances, stable rotational wave packets and strong sen- 
sitivity of their coherent superpositions to the phase re- 
laxation width. Note that quasiperiodic structures in the 
excitation functions are present also for other heavy-ion 
colliding systems . We have also described thoroughly 
how to extract the time power spectrum from the data 
on the excitation function oscillations. 

It should be mentioned that a conceptually different 
optical model description of the 12 C+ 24 Mg elàstic scat- 
tering at backward angles was presented in Ref. [Í3| . 
where the quasimolecular structures were ruled out. This 
conclusion was based on a qualitative description of an- 
gular oscillations at backward angles for only three ener- 
gies within the same energy range which is the subject 
of our study in this paper. Unfortunately, that analysis 
is inconsistent with the presence of fine energy structure 
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of ~ 0.3 MeV at = 180° ÜJ, which unambiguously re- 
veals the presence of time-delayed processes at backward 
angles. On the contrary, the optical model description 
relies on the energy average 5-matrix and therefore is 
cxpected to reproduce broad ~ 2 — 3 MeV energy struc- 
tures in the cross section. From our point of view, the 
fact that the presence of the fine energy structures was 
not taken into account led to an incorrect qualitative and 
then quantitative identification of the physical picture of 
the collision process in Ref. [Í3| . 

In conclusion, we have demonstrated a strong sensi- 
tivity of Schrodinger cat states to the phase-relaxation 
width in complex quantum collisions. This should per- 
mit to determine the phase relaxation time from mea- 
surements of the excitation functions for 12 C+ 24 Mg elàs- 
tic scattering at backward angles. Such an experiment 
would be desirable since the Schrodinger cat states pre- 
dicted in [ïg| for 12 C+ 24 Mg scattering involve - 10 a -10 4 



many-body configurations of the highly excited interme- 
diate complex. To the best of our knowledge, the inter- 
nal interactive complexity of these quantum macroscopic 
superpositions dramatically exceeds |19| all those previ- 
ously experimentally realizcd. The proposed method can 
also be applied for determination of the phase relaxation 
time from the data on excitation functions for bimolecu- 
lar chemical reactions. 
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